NASA TM X-52291 


NASA TECHN ICAL 

MEMORANDUM 


NASA TM X-52291 



THE SELECTION OF SIMPLE MODEL ATOMS 

FOR CALCULATIONS OF ELECTRON DENSITY 

IN NONEQUILIBRIUM, LOW TEMPERATURE Cs PLASMAS 


by John V. Dugan Jr. 
Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PAPER proposed for presentation at 27th Annual 
Conference for Physical Electronics sponsored by 
Massachusetts Institute of Technology 
Cambridge, Massachusetts, March 20-22, 1967 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D.C. • 1967 


THE SELECTION OF SIMPLE MODEL ATOMS FOR CALCULATIONS OF 


ELECTRON DENSITY IN NONEQUILIBRIUM, LOW 

TEMPERATURE Cs PLASMAS 

by John V. Dugan Jr. 

Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PAPER proposed for presentation at 

27th Annual Conference for Physical Electronics 
sponsored by Massachusetts Institute of Technology 
Cambridge, Massachusetts, March 20-22, 1967 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


E-3904 


THE SELECTION OF SIMPLE MODEL ATOMS FOR CALCULATIONS OF ELECTRON 
DENSITY IN NONEQUILIBRIUM, LOW TEMPERATURE Cs PLASMAS 

John V. Dugan Jr. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Free electron number densities N e are calculated for cesium plasmas with low electron tem- 
peratures T e (0.15 - 0. 3 eV) using 5-, 3-, and 2-level atomic models. These Ne values are ob- 
tained for a Maxwellian distribution of free electrons; this distribution was shown to be a good 
approximation for T e > 2100° K in optically thick Cs plasmas. The sensitivity of N e to the 
choice of electronic levels, statistical degeneracy and radiative capture coefficients is presented 
for plasmas either optically thin or thick to resonance radiation. The steady state equation for the 
electron number density is discussed in the limits of low and high neutral density N^ s . The ion- 
ization fraction f - N 6 /Nq s was studied as a function of Nq s for the range of T e from 
1740° K {0.15 eV) to 3000° K. Results are presented as a function of T e for Nq s = 3x10^ m“®, 
(typical values for certain thermionic diodes). Both sets of results are compared with the cor- 
responding Saha curves for plasmas in local thermodynamic equilibrium (LTE). A dis- 

cussion of rules for choosing reasonably accurate but simple atomic models is included. These 
electronic levels which serve as free electron source terms are identified as functions of N^j s 
for fixed T e . Even a 5-level model gives N e values which are only in fair agreement with, 
optically thick results for a 26-level model under certain conditions. The lumped level degen- 
eracy gj^ can be made a function of N^, g and T e to bring the M = 5 results into better agree- 
ment with the more reliable M = 26 results of Norcross and Stone. However, in the absence of 
such a comparison it does not appear that construction of accurate but simple atomic models is a 
straightforward process. 


INTRODUCTION 

There is much interest in the accurate calculation of free electron number density N e in non- 
equilibrium plasmas of low electron temperature T e that exist in Cs plasma diodes. The pur- 
pose of this study is to arrive at rules for constructing simple but satisfactory model atoms. Such 
models must ensure reasonably accurate values of N e for the aforementioned plasmas as well as 
allow simplified calculations of the free electron distribution function. The accuracy of such 
models depends upon the number and arrangement of electronic energy levels as well as inelastic 
cross sections, radiative lifetimes and radiative capture coefficients. For the steady state plas- 
mas of interest, the free electrons are assumed to have a Maxwellian distribution. Solutions of 
the Boltzmann equation (including all elastic and inelastic collision terms) along the lines described 
in Ref. 1 indicate that the distribution function f e (u) is Maxwellian down to T e values of 2100° K 
(for Cs plasmas optically thick to resonance radiation). Although the numerical method described 
in Ref. 1 was convergent only for T e s 2300° K, subsequent calculations down to T e = 2100° K 
have been found to yield essentially Maxwellian f e (u). The high energy tail of f e (u) may be non- 
Maxwellian for T e values <2100° K but the N e values can be relatively insensitive to this tail 
{u > 3. 89 eV) for optically thick plasmas. This insensitivity is present where ground state ion- 
ization becomes a relatively unimportant rate process. 

DETAILS OF THE CALCULATION 



A model atom consisting of M levels is studied. The rate of change of the number of bound 
electrons in any level L is^: 
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In Eq. (1) Kl_*k (m** sec -1 ) is the excitation coefficient for the collision induced transition L-»K 
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Kk-»l (m a sec - -*-) is the corresponding coefficient for de-excitation. The quantities 
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(m° sec”-*-) is the ionization coefficient for level L. (sec -1 ) is any 

radiative transition probability for de-excitation into the Lth state and Aj^ is the sum of the 
Al_*k (sec -1 ) for radiative transition to all levels below L. 

For the steady state treated here Nl = 0 and the number densities must satisfy the plasma 
normalization condition: 
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Equations (1) and (2) together with the condition of charge neutrality (N e = Nj), serve to determine 
N e and the Nl values for a specified f e (u) and initial Cs number density Nq s . Equation (2) can 
be rewritten as: 

M 

f + N L = 1 (2a) 

L=1 

where f is the ionization fraction and is the normalized (to Nq s ) population of the Lth state. 
In the steady state the total ionization rate from bound levels must balance the total capture rate of 
free electrons into all such levels. This relation provides a useful check on the numerical calcu- 
lations. Its application indicated that accurate numerical solutions for f were not possible at low 
values of N e without tight convergence criteria. 

The condition, N e = 0 can be written as: 

M M 

[(K° a P)*f 2 + jS L fj = ^ K^ n (3) 

where (K ca P)* (m 3 sec -1 ) is an effective two-body capture coefficient = N^gK^ 3 ’^. Thus the 
equation for the ionization fraction is essentially quadratic; however, the normalized populations 
Nl are strong functions of N^ s and optical thickness for fixed T e . The ionization coefficient 
K L, n — <Q ion (v e )v e ) is of course a function of T e alone. 

MODEL ATOM PROPERTIES 

The optimization of the model atom depends upon the choices of cross sections for electron im- 
pact and coefficients and frequencies for radiative processes involving the M 1 * 1 level. This level 
must be adjusted so as to simulate the presence of many missing excited levels 1 . The assign- 
ments of ionization potential E m statistical degeneracy g M and Aj^ and /3 M values thus be- 
come ciiiicai for accurate calculations of f values. The N]y[ for this lumped ievei is a com- 
plicated function of the level structure for variable N^ s (see Eq. (1)). The schematic diagram 
for the three models is shown in Fig. 1 with values of El, Al-»-k a* 10 * SL‘ 

The early form of the Gryzinski cross sections has been used to compute ionization coefficients 
and excitation coefficients for optically allowed transitions only. 2 The Gryzinski exchange formu- 
lation was used to compute excitation coefficients for optically forbidden transitions. 3 Although 



the exchange cross sections have relatively small values at maximum, their slopes are steep just 3 
above threshold. Since the excitation and ionization thresholds are usually »kT e , the behavior of 
the cross section in this region determines the value of the excitation coefficient. The lumped 
level was treated as an optically allowed state for all excitation (and de -excitation) collisions. 

Most allowed and exchange cross sections (monoenergetic) give reasonable agreement with ex- 
periment over a range of electron energies from threshold to maximum^. The experimental value 
of the most important excitation cross section ^ as a s l°P e roughly 3 times the Gryzinski 

value at threshold^. 

The radiative capture coefficients 0 l were taken from Ref. 5 where they were calculated 
using an adjusted quantum defect method. These calculations agree well with known oscillator 
strengths and recombination cross sections. 

The radiative transition probabilities were taken from Ref. 6; the agreement of those calcula- 
tions with cited experimental values is within 50% for all important values (see Ref. 7). 

The most critical Alk value is for the 2-*T resonance transition; this line is strongly absorbed 
for all experimental plasmas of interest®. For M = 3 the A32 value was set equal to the sum of 
A52, A42, and A32. The value of A21 was used for all three optically thin cases. 

RESULTS - OPTICALLY THIN PLASMAS 

Since the only available M = 26 results (optically thin) are several points for T e ~ 3000° K, 
the M = 2 and 3 results will be compared with 5-level results for those plasmas. Since ground 
state ionization is a significant process at low Nq s values, however, a non-Maxwellian tail in 
fg(u) could cause large changes in calculated values of f in this region. 

Low N^g limit. - In the limit of low Nq s values three -body recombination is negligible and 

the free electron number density is determined by (from Eq. (3)) 
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In the optically thin case, radiative capture into excited states is the important capture process 
which balances ground state ionization. Then f is given approximately by: 


f « K 


ion 
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which simplifies for M = 2. The product Kj° n I r j_= Kj° n (for allM), however, since the ground 
state population is greater than 0.95 N^ s for optically thin plasmas in the low N^ s regime. The 
value of f is then simply determined by the ratio of the ground state ionization coefficient to the 

sum of the largest radiative capture coefficients. For cesium the 0g (5D level) value dominates 
M 

the ^ /3 l term; it is three times the maximum 0 value for other excited levels, 0(6P), and 
200 times the 0^ (6S) value 5 . Thus, in choosing simple model atoms, it is important to assign 
a large value of 0jy{ to simulate the missing excited levels at low Nq s . The 0 5 assigned for 
M - 5 is the sum of /3 values for the six excited states: 7P, 6D, 8S, 4F, 8P, 7D . 

The results of f versus N^g for 2-, 3-, and 5~level models are shown in Fig. 2 for an elec- 
tron temperature of 3000° E (0.26 eV). The f values for M = 2 (6S and 6P levels) with 0g = 
0(6P) + 0(5D) agree with the M = 5 results within 50% at low values of N^ s (^3xl0^° m~3). 

This agreement can be improved by adding the sum (/3(7S) + /3 5 ) to the 02 value. For a given 
excitation cross section Ql-^k t ^ ie de-excitation coefficient Kj^-l is inversely proportional to 
the upper state degeneracy gj£. The L^ 1 level 3 -body capture coefficient is directly propor- 
tional to the g L value. The change of degeneracy for the M^ 1 state (M = 2, 3) does not affect 
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the f value at low Nq s because the value of is unchanged. However, the degeneracy does 
affect the maximum ionization fraction calculated for the M = 2 and 3 models. This latter be - 
havior can be described in terms of the free electron source terms although there is no simple 
limiting expression. The degree of departure of f from LTE conditions is evident by direct 
comparison with the Saha curve in Fig. 2. 

High Nc s limit. - In the limit of high N^g values 3~body capture becomes the dominant loss 
mechanism for free electrons. Radiative capture is negligible and the value of f can be approxi- 
mated by: 



The excited levels have the largest ionization coefficients since the ionization cross section Ql * 1 
is a strong inverse function of ionization potential El. These levels are most nearly in equilib- 
rium with the free electrons and f approaches the Saha result. For the simplest model atom, 

M = 2, one must choose values of gjy[ and E^ so as to preserve the "correct" ratio of ioniza- 
tion rate to 3-body capture coefficient. In this study the value of E2 was fixed at the 6P value, 
2.43eV, and the value of g2 was varied for best agreement with M = 5 results. However, it 
will be shown that the Mth state for M = 5 is not so simply adjusted to match the M - 26 results. 

For a given ratio of ionization and capture coefficients for the lumped level N2 must vary with 
Nq s so as to account for the populations of the missing excited levels near the Saha limit. For 
M = 2, N2 increases so rapidly (oc N^ s °- 95 ) that K| on N2 exceeds (K| a P) + f 2 at intermediate 
values of Nq s . The f is relatively small and only a gradual function of N^s so the 3-body rate 
does not become important for the optically thin case until high Nq s values. The maxima in the 
f curves of Fig. 2 occur approximately where the collisional lifetime = ^2^e)~^ equals the radi- 
ative lifetime (A21) -1 for the first excited state (see Ref. 9). 

The M = 3 and 5 cases are similar in that net ionization from excited states is a significant 
process for Nq s > 10 22 m~ 3 . The M = 2 model is unique because the ground state always serves 
as the free electron source term (i.e. exhibits net ionization). Thus in this case the lumped level 
necessarily becomes a net capturer since N e -- 0. The large /3 2 value (60 times ( 3 -^) ensures 
that this level remain a net capturer at low values of Nq s . The sources of ionization for the 
M = 3 and 5 models can be determined from the net ionization rates. The rate for the L'h state 
can be written from Eq. (3) as: 

R net = K g n Nl - [(KffP)*f 2 + 0 L f] (6) 

Only R^ et values which are positive will be considered, i. e. , states which are free electron 
source terms. The variation of such source terms (normalized to the sum of positive terms) 

Rl is shown for T e = 3000° K, optically thin, in Fig. 3 for M = 5 and 3. The L = 4 level 
(7S) makes the largest contribution to the free electron population at T e = 3000° K for M = 5, 
Nq s > 2X10 2 * m -3 . This behavior differs slightly from the T e = 0. 2 eV results where the L - 3 
state (5D) has the largest Rl value at high Nq s . 

The results for M = 5, g 5 = 50, E 5 = 0. 6 eV of Fig. 2 are in good agreement (within a factor 
of two) with Ng/^NgJgaha results for the 26-level model at T e = 3000° (see Table I). No 
optically thin results for M - 26 were available for comparison at lower T e values. The sim- 
pler models should give good agreement with the M = 26 case at low (— 10-*- 9 m -3 } since the 

lumped /Ijyj values have been chosen so as to produce that behavior 

RESULTS - OPTICALLY THICK PLASMAS 

Low Nq s limit. - At low values of N^s ( <10 19 m~ 3 in Fig. 4(a)) f is determined by the bal- 
ance between ionization from excited states and radiative capture. Just as in the optically thin 
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case for high Nr-,, the excited states are free electron source terms, but at low Np s . The 
dominant is R$> for Nq s — 10 iy m”* 1 in the 3- and 5-level models. Only the ground state 

has a positive Rj^ value for M - 2 through the range of Nq s values. The approximate rela- 
tion for the ionization fraction becomes: 







For M - 5 the above ratio should give reasonably accurate f values when just the L - 2 and 3 
ionization terms and the /3g value are included. However, at intermediate values of N^ s 3-body 
capture is comparable with the radiative process and the lumped level must represent the missing 
excited states. The excited states are free electron source terms at low N^ s because their 
populations (especially N 2 ) are enhanced in the absence of radiative de-excitation. The ioniza- 
tion rates become large even when populations Nl, L > 1 remain « 10 " 2 -*• 10"® of Nj because 
the ionization coefficients K^ n » K^ on . The monoenergetic cross sections Q^ n , L > 1, have 
lower energy’ thresholds and much higher slopes than Qj° n in the threshold region. The excitation 
coefficients for collision-induced transitions between excited levels are correspondingly larger 
than the ionization coefficients and have low thresholds (some <1 eV). 

The f values for the optically thick case are plotted versus Nq s for T e -- 3000° K in 
Fig. 4(a). The lumped value should guarantee accuracy of the simpler models at low n£ s . 
This is true if the populations N 2 have nearly the same values as for M =■ 26. Decreasing /3 m 
for M = 2 gives better agreement with M = 26 at N^ g = 10 2 ® m ~ 3 but the asymptote is incor- 
rect at 10 1S m“ 3 . For M = 3, the f value is increased at all Nq s by making the plasma thick 
to all line radiation and lowering g^j from 50 to 10. There is no peak in f versus Nq s for 
M - 2, g 2 = 50 since the L - 2 state has a large 3-body capture rate at intermediate N^ s due 
to its large g^ value and relatively high N e . Lowering the degeneracy to g 2 = 6 (not shown) 
does not change the curve shape since 3-body capture continues to dominate excited state ioniza- 
tion. Decreasing the gjyj from 50 (curve 4) to 10 (not shown) for M - 3 does produce a maxi- 
mum in f because 3-*- 2 de -excitation is increased enough to raise the N 2 value significantly. 
This L = 2 state is the main free electron source term for M - 3. 

The f values for the optically thick case, T e = 2321° K, are plotted versus in Fig. 4(b). 

Just as for T e = 3000° K the f is determined by ionization from the ground level for M ~ 2 
and from excited states for M - 3 and 5. These ionization processes balance radiative capture 
into all levels and the curves remain flat where the populations of the ionization sources remain 
constant, just as for 3000° K. The features of the curves are identical to the T e - 3000° K re- 
sults throughout the range of N^g values. The peak for M-5 is again caused by ionization 
from the higher excited states L = 3 and 4. Four points of M = 26 results are also shown in 
Fig. 4(b). Comparison with the latter results indicate that the simpler atomic models overesti- 
mate departures from LTE at values of Nq s below 3xl0 22 m -3 . It should be noted, however, 
that the N e values from M - 26 are considerably lower than Saha values for N^ s ^ 10 2 * m' 3 . 
The ionization fractions are approximately 10 ” 2 of the 3000° K values because excited state 
populations and ionization coefficients are considerably lower. The limiting (low N^ g ) value of 
f for M = 3 is only 10 " 2 of the f computed for a 3 -level model in Ref. 9. For M = 2 the value 
of f is about 1/15 of that for a 2-level model in Ref. 9. These large differences are due to the 
present use of more reliable excitation cross sections and and values. 

The variation of the free electron source terms for M •- 5 at T e = 2321° K (0. 2 eV) and 
3000 K are shown in Figs . 5 (a) and 5(b). respectively Qualitatively, the dependences of 
are similar but the region of rapidly changing contribution is shifted by a factor of 10 2 to higher 
N 0 g values at the lower value of T e The M - 26 results of Ref. 7 indicate that the maximum f 
is nearly 5 times the M = 5 value for Nq s = 10 21 m~ 3 , T e - 2321° K. Figure 5(b) shows that 
the main contributor of free electrons is the L - 3 state for N^ s > 3x1 0 2 ^ m ' 3 and the ground 
and first excited states for Nq s values <3xl0 2 ^ m~ 3 . If additional excited states are to become 
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new ionization source terms, their relative populations must fall off slowly enough with binding 
energy that their higher ionization coefficients can raise f . The sudden drop in the L = 3 contri- 
bution in the region Nq s £ 10 21 m -3 is a result of the large value since radiative capture 
becomes the dominant electron loss process in that region (N e ^ 5X10-*- 3 m -3 ). 

The agreement between the M = 5 and 26 results improve at high N^g (= 10 22 m -3 ) for 
T e = 3000° K as can be seen in Fig. 4(a). The f for M = 5 is about 90% of the f value for 
M = 26. Also, the disparity near peak (Nq s = 10 2 ® m -3 ) is explainable in terms of the value 
just as the T e = 2321° K results. Figure 5(b) shows that the L = 4 level is the chief contributor 
of free electrons for N^ s > 3x10^ m -3 . This occurs because the "ladder” ionization mechanism 
operates if the excited state populations are large enough. For optically thick plasmas this is the 
case. It is plausible, therefore, that including additional excited states would raise the value of f 
in this T e region if these states serve as large source terms with Rl values = R^. 

High Nq s limit. - At high values of Nq s the excited levels are still the free electron source 

terms. Again the balance which determines f is between excited state ionization and three-body 
capture into these same levels. The f values, optically thick, are uniformly larger at lower 
Nq s values than in the optically thin case. This is because the important Nl approach their 
Saha values at lower N^s (— 3xl0 2 * m -3 ). The value of f can be satisfactorily calculated from 
Eq. (5) including terms, L = 2 to M. 

Results for T e = 1740° K (0. 15 eV). - Curves of f versus N^« s were also calculated at 
T e - 1740° K. There is no structure in these curves for N^, g < 10 2 ^ m -3 , optically thin and 
Nq s < 10 23 m -3 , optically thick. One point was compared with the M = 26 results at N^g = 
lO 2 ^ m~ 3 , optically thick. The M = 26 model gives an f which is 3 times the value for M = 5, 
g5 = 50. Since the chief free electron source term is L = 5 for M = 5 the difference between the 
M = 5 and 26 results becomes very sensitive to the lumped state. This level must now simulate 
the source terms for N e . 

The simpler models do not agree well with each other at T e = 0.15 eV, optically thick. This 
behavior is shown in Fig. 6 where f is plotted as a function of T e at n£- s = 3X10 23 m -3 for 
M = 2, 3, and 5. The largest disparity occurs in the optically thick case where the f for M = 2 
is only 2% of the value calculated for M = 5. This result is consistent with the fact that the 
former model does not allow for excited state source terms. The M = 3 model gives an f value 
midway between the M = 2 and 5 results. The values of f, optically thin, are only = 10 ~ 3 of the 
optically thick values; this is because the ground state is the sole ionization term in the former 
case. Thus at this low T e value the ionization mechanism changes from simple ground state ion- 
ization to a "stepwise” ladder process (for N^j g = lO^ 3- *’!© 2 ® m -3 ) as the optical thickness 
changes. As the electron temperature is increased this disparity between optically thin and thick 
results decreases markedly. For fixed values of /3]y[ and Ej^, agreement with M = 26 results 
for simpler models can be improved by using variable g-^ = F(Nq s , T e ) or making the plasma 
thick to all radiation. However, any scaling laws for construction of atomic models remain to be 
investigated in the absence of more detailed results for comparison. 


CONCLUDING REMARKS 

For certain values of electron temperature and neutral density 2- and 3-level models can be ad- 
justed to give N e values which agree with results for 5-levels. However, for the experimentally 
interesting case of optically thick Cs plasmas, results for a 26 level model^ indicate the "best" 

M = 5 model can seriously underestimate the ionization fraction t. This is especially true fur 
values of electron temperature <3000° K. This disparity occurs because the sources of free elec- 
trons are highly excited levels as a result of a "ladder" ionization mechanism. It does not appear 
likely that convenient adjustment of simpler models will improve agreement with the M = 26 re- 
sults. For proper choices of Ejy[, c^l for the lumped level a variable gM(Ncs) will pro- 
vide satisfactory agreement for fixed T e . However, it appears that increasing the number of 



levels is a more desirable approach than this to ensure calculation of accurate Ne values. Fur- 
ther results for M = 26 are needed to evaluate the simpler model atoms for the low temperature 
(1740° - 2321° K) optically thin plasmas. 


TABLE I. - (T e = 3000° K; OPTICALLY THIN) 


M 

Ne/(N e ) S aha 

5 

2.4X10- 2 

6.8x10-2 

0.35 

0.75 

26* (Ref. 7) 

=1.3xl0- 2 

=6. 5x10 “ 2 

= . 32 

= .75 


♦Points estimated from plots without grid. 
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Figure 3. - Plots of normalized net ionization 
rates as a function of Cs neutral density in 
plasmas optically thin to resonance radiation 
for 3- and 5- level atomic models. T e ■ 3000° K. 



Figure 2. - Plots of ionization fraction as a func- 
tion of Cs neutral density in plasmas optically 
thin tc resonance radiation for 2 ~, ^ nri 

5- level atomic models. T e = 3000° K. 



(a) T e ■ 3000° K. Two points for a 26- level 
model are shown for comparison. 


Figure 4. - Plots of ionization traction as a func- 
tion of Cs neutral density in plasmas optically 
thick to resonance radiation for 2 ~, 3-, and 
3- level atomic models. 
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Figure 4. - Concluded. 



(a) T e = 3000° K. 

Figure 5. - Plots of normalized net ionization 
rates as a function of C:> iieuti'a' density in 
plasmas optically thick to resonance radiation 
for the 5- level atomic model. 




Figure 6. - Plots of ionization fraction as 
a function of electron temperature in 
plasmas optically thin and thick to reso- 
nance radiation for M = Z~, 3- ( ana 
5- level atomic models. 



